Solid phases and pairing in a mixture of polar molecules and atoms 
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We consider a mixture of hard core bosonic polar molecules, interacting via repulsive dipole-dipole 
interaction, and one atomic bosonic species. The mixture is confined on a two-dimensional square 
lattice and, at low enough temperatures, can be described by the two-component Bose-Hubbard 
model. The latter displays a extremely rich phase diagram including solid, superfluid, supersolid 
phases. Here we mainly focus on the checkerboard molecular solid, stabilized by the long range 
dipolar interaction, and study how the presence of atoms affects its robustness both at zero and 
finite temperature. We find that, due to atom- molecule interaction, solid phases can be stabilized at 
both, (much) lower strengths of dipolar interaction and higher temperatures, than when no atoms 
are present. As a byproduct, atoms also order in a solid phase with same melting temperatures 
as for molecules. Finally, we find that for large enough interaction between atoms and molecules a 
paired supersolid phase can be stabilized. 



Recent years have witnessed spectacular advances in 
atomic and molecular physics. Ultracold bosonic atoms 
have been successfully loaded in one- 0, two- three- 
dimensional Q optical lattices with full control over the 
hamiltonian parameters. The Mott-insulator (MI) to 
superfluid (SF) transition has been experimentally ob- 
served and extensively studied theoretically [J, [f| ■ More- 
over, more recently, new experimental efforts towards 
trapping two component atomic mixtures in optical lat- 
tices have been put forward [fj|-flo| . In the presence of a 
second component, new fascinating phenomena and ex- 
otic quantum phases which cannot be accessed with sin- 
gle species gases, become available (Til 

Ultracold polar molecules, interacting via long range 
dipolar interaction, also represent a fertile ground for 
the study of novel exotic phases of matter. The dipole- 
dipole interaction can be tuned and shaped via external 
static and microwave electric fields [l3j]. Self-assembled 
dipolar crystal, dipolar Bosc-Einstein condensate, super- 
solids, are among the phases which can be stabilized with 
cold polar molecules (Til [l5j. Besides, interesting appli- 
cations in fundamental physics, e.g. measurement of per- 
manent electron dipole moment, in quantum chemistry 
and few-body physics, and in quantum information pro- 
cessing have been envisioned [161 ]. 

Observation of many body quantum phases requires 
high phase space density of molecules in stable ground 
states. Cooling high density samples of molecules to 
quantum degeneracy has proved to be challenging due to 
their complex internal state structure. Only very recently 
samples of high density polar and non-polar molecules 
in the rovibronic ground-state and a selected hyperfme 
state were realized [Tt], [H| . The Jila group has also re- 
ported on observations of dipolar collisions and evidence 
of spatial anisotrpy of dipolar interaction in a sample of 
KRb molecules with non zero induced dipole moment, 
i.e. aligned by external electric field (l9| . The Inns- 
bruck group has produced a high density sample of Cs2 
molecules in their rovibronic ground state trapped into 
an optical lattice [TU . The procedure used in [l8| can be 
generalized to heteronuclear molecules. 



In the present work we consider a novel system consti- 
tuing of a mixture of polar molecules and one atomic, non 
reactive species. This system is realizable experimentally 
since residual atoms naturally occur after molecule for- 
mation is completed and can be selectively removed [2(| ■ 
We are interested in studying how the presence of the 
atomic species affects molecular solid phases stabilized 
by dipolar interaction [l5[. Hence we consider the case 
where only reactive atoms are removed from the system, 
and both, molecules and atoms, are in the absolute (in- 
cluding hyperfme) ground state, so that inleastic colli- 
sions are forbidden. Examples of non reactive mixtures 
include (KRb+Rb), (LiCs+Cs), (RbCs+Cs) 

We consider bosonic polar molecules in a static elec- 
tric field along the z-direction which induces a dipole 
moment d parallel to the field. Molecules interact via 
dipole-dipole interaction V dd (R) = . If the 

system is confined on the xy plane (by e.g. a strong 
optical lattice potential along the z-drection) in order to 
avoid instabilities which might arise due to the attractive 
part of the dipolar interaction, the resulting effective 2D 
potential is isotropic, ~ fj. An additional 2D opti- 
cal lattice with spacing a (which we use as unit length) 
further confines the system. If atoms are present, they 
interact with molecules via Van dcr Waals (anisotropic) 
interaction, i.e. Vydw( r ) = %, which we describe via an 
effective s-wave scattering length, and via dipole-dipole 
interaction V^'L mol (r) = ^f-. Since atoms have small 
polarizabilitics, dipolar interaction between atoms and 
molecules can be neglected, provided the electric field 
is small. Indeed, for DC electric fields Edc ~ kV/cm, 
and molecular (induced) dipole moment d ~ ID, the 
dipolar interaction at a typical lattice spacing distance 
a ~ 500 nm is V£j~_ mol (a) ~ 0.1 — 1 Hz h (to be compared 
with molecule-molecule dipolar interaction V dd (a) ~ kHz 
h). At low enough temperature, the system is described 
by the two-component Bose-Hubbard Hamiltonian: 

H=-J b J2 4b, Ja £ aja, + v£ ^ 
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-U ab J2n> b i + {£u aa - 1) 

where &i (oj) and b\ (af) are bosonic operators for 
molecules (atoms) with b\ 2 = (the hardcore limit is well 
justified if one starts with low enough molecular density 
so that there are no doubly occupied sites [HI), and n*' b 
is the occupation number operator; the first and second 
terms in Eq. Q] describe the standard kinetic energy with 
hopping rate J a ,b for atoms and molecules respectively, 
the third term describes the repulsive dipole-dipole in- 
teraction between molecules with = \i — j\, the fourth 
and fifth terms describe the atom-molecule and atom- 
atom contact interaction, respectively; fi a ,b is the chemi- 
cal potential. 

We study model (JlJ via large scale Monte-Carlo sim- 
ulations by the Worm Algorithm [23[. An Ewald 
summation is used in the interaction potential so that no 
cutoff in the range of the dipole-dipole interaction is in- 
troduced. The complete phase diagram is extremely rich. 
It includes mixtures of a molecular solid or supersolid 
(SS) with MI or SF of the atomic component, double SF 
phases, etc. Besides, one expects the presence of paired 
phases which have been proved to exist in multiconrpo- 
nent atomic systems (see e.g. [IJ E3)- I n the present 
work we mainly focus on the checkerboard (CB) molec- 
ular solid phase at density rib = 0.5 0]. Our results can 
be summarized as follows. Due to atom-molecule inter- 
action, the CB molecular solid appears at weaker dipo- 
lar interaction compared to the purely molecular sample. 
Correspondingly, solid melting temperatures are higher 
than in the absence of atoms. As a byproduct, atoms are 
also ordered in a CB solid phase with same melting tem- 
peratures as for molecules. We also find that, for large 
enough interaction between atoms and molecules a paired 
supersolid (PSS) phase, i.e. a SS phase of a composite 
(atom+molecule) object, can be stabilized upon doping 
of the solid. 

Let us first discuss the parameter space we focus on. In 
the following we consider J a = Jb = J (although the re- 
sults discussed below hold in general for J a ~ Jb). This 
can be easily achieved by making use of two indepen- 
dent lattices for atoms and molecules. We work at den- 
sities n a ~ rib ~ 0.5, in the limit of hard-core molecules 
and atoms, and consider attractive interaction between 
atoms and molecules, i.e. U a b > 0. Results still hold for 
repulsive atom- molecule interaction, i.e. U a b < 0, pro- 
vided one replaces atoms with holes (and viceversa) in 
the following discussion. The hard-core limit is naturally 
reached for molecules [22[, as we have mentioned above. 
For atoms, it can easily be reached via Feshbach reso- 
nances. The qualitative results found here, though, are 
still valid for soft core atoms as long as U aa > \U a b\ J- 
For example, at U aa ~ 2.5U ab , and U a b = 8J no qual- 
itative difference from the hard-core case was observed. 
This can be understood from the fact that, in the limit 
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FIG. 1. (Color online) Checkerboard lobe with parameters 
n a = 0.5, J a = Jb, and U a b/J = 5, 8, down traingles and 
diamonds respectively. The lobe in absence of atoms is also 
plotted for reference (circles), with data taken from [Q. Er- 
rorbars are within symbol size when not visible. 
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FIG. 2. (Color online) Checkerboard lobe with parameters 
n a = 0.5, J a — Jb, and U a b/J = 12, 20, up traingles and 
squares respectively. The lobe in absence of atoms is also 
plotted for reference (circles), with data taken from [l5|]. Er- 
rorbars are within symbol size when not visible. 



Uaa ^ \U a b\ ^ J, the difference between hard- and soft- 
core appears in fourth order perturbation theory. The 
constraint on particle numbers can be satisfied in the 
presence of an external harmonic confinement, which pro- 
vides a scan over chemical potential. We have indeed ob- 
served that, when molecules form a solid, the atomic den- 
sity follows the molecular one for a wide range of atomic 
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chemical potentials and already at relatively low atom- 
molecule attraction, implying that n a ~ n& ~ 0.5 will be 
easily satished in some regions of the trap. 

Let us now turn to the results. In Fig. [1] and [5] we 
study how the onset of the CB solid is affected by the 
presence of atoms, upon varying the interaction between 
atoms and molecules limiting ourselves to values U a b / J 
for which atomic density also orders, forced by the molec- 
ular solid. In Fig. [1] we plot the CB lobes for U a b/J = 5 
(down triangles), U a b/J = 8 (diamonds), and, for refer- 
ence, in absence of atoms (circles). (We do not investi- 
gate the nature of phases outside the lobe, as establish- 
ing the full phase diagram of model ([T]) is beyond the 
scope of this work). The boundaries of the lobes are de- 
termined by monitoring n&(/Xb) for system sizes L=12, 
16, 20, 24, 30 (with extrapolation to infinite size when 
necessary). With no atoms in the system the solid is 
stabilized for dipolar interaction strength V/J > 3.6. 
In the presence of atoms, instead, the solid already ex- 
ists for V/J > 2.5 (1-75) at atom- molecule interaction 
Uab/J = 5 (8). Similarly, Fig. [2 shows CB lobes for 
Uab/J = 12 (squares), U a b/J = 20 (up triangles), and in 
absence of atoms (circles). Here the effect of the atom- 
molecule interaction is even more dramatic. The solid is 
stabilized for V/J> 0.5, a factor of 7 smaller than the 
case of no atoms present in the system! The renormal- 
ization of the tip of the lobe towards smaller values of 
V/J can be understood by noticing that, for U a b/J val- 
ues considered here, atoms also ordered in a CB solid. 
One can therefore visualize a solid of a composite object, 
constituent of one atom and one molecule, which has an 
effective hopping J e s oc -g— < J responsible of rescaling 
the onset of CB solid towards smaller dipolar interaction 
strengths. 

The results summarized in Fig.s [1] and [2] have impor- 
tant consequences for experiments, as they show that in 
presence of non reactive atoms, smaller induced dipolc 
moment is needed in order to stabilize the solid. If 
one considers (RbK+Rb) mixtures loaded in lattices of 
depth V a ~ 8E? and V b ~ 7E b r (E™ is the recoil en- 
ergy of species a) and with spacing a ~ 500nm, one has 
J a ~ Jb ~ 80Hz h and U a b/J ~ 8, under the assump- 
tion _ a KRb-Rb ^ I50a - In absence of atoms, observa- 
tion of CB solid requires dipolar interaction V > 3.6 J ~ 
288Hz h which can barely be reached with fully polirized 
molecules. In the presence of atoms, instead, one needs 
V > 1.5J <~ 120Hz h, which requires an induced dipolc 
moment d > 0.351? (easily achievable in the laboratory 
with alkali dimers). 

The enhanced robustness of the CB solid further man- 
ifests itself in larger melting temperatures compared to 
the case of no atoms present. In Fig. [3] we plot the criti- 
cal temperature for the melting of the CB solid T C J J as 
a function of V/J for U a b/J = 5 (diamonds), 8 (down 
triangles), 12 (squares), 20 (up triangles), and with no 
atoms (circles). We find that critical temperatures have 
already saturated at U a b / J = 8 and can be up to a fac- 
tor of 2 larger than in absence of atoms. (We have not 




FIG. 3. (Color online) Checkerboard solid melting tempera- 
tures with parameters n a — nt — 0.5, U a b/J = 5, 8, 12, 20 
diamonds, down triangles, up triangles, squares respectively 
and, for reference, in absence of atoms (cricles). 

resolved T c /J close to the tip of the lobes where sim- 
ulations are more challenging.) Within our statistical 
errors, melting temperatures for the atomic density wave 
are the same as for molecules, suggesting that a solid 
of a composite (atom+molecule) object survives up to 
the melting transition line. To appreciate the effect of 
atoms let's consider an induced dipole moment of 0.5 D 
in a (RbCs+Cs) mixture loaded in a lattice with spac- 
ing a = 400 nm and depths V a{b ) ~ 9(7)E? [b) . With 
— a s ~ 150ao we obtain J a ~ Jb ~ 70Hz h, U a b/J ~ 12, 
and a solid melting temperature 2 times larger than in 
absence of Cs atoms. 

Next we turn to the discussion of paired phases. Here, 
we would like to focus on the possibility of realizing a 
paired supersolid (PSS), i.e. a SS of a composite ob- 
ject formed by an atom and a polar molecule (we remind 
the reader that a SS is a state of matter featuring both 
supcrfluid and crystalline order). In order to check for 
the existence of a PSS we have calculated the structure 
factor and the winding number statistics, upon doping 
the CB solid with interstials, i.e. we work at densities 
n-a ~ Tib < 0.5. The structure factor S(k), which charac- 
terize solid phases, is defined as: 

S a (k) = ^cxp[ik(r-r')](n Q (rK(r')>/iV Q (2) 
r,r' 

with a = a, b, k = (%, n) the reciprocal lattice vector for 
the CB solid, and N a the total number of particles of 
species a. Paired superfluidity (PSF), i.e. superfluidity 
of the composite object, is instead identified by simul- 
taneuos finite statistics of the sum of winding numbers 
squared < (W a + Wb) 2 > oc n PSF , where n PSF is the 
supcrfluid stiffness, and zero statistics of the difference 
of winding numbers squared < (W a — Wb) 2 >■ We 
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FIG. 4. (Color online) Structure factor Sb(ir,Tr) (right axis) 
and pair superfluid stiffness n^ SF (left axis) as a function of rib 
with parameters J a — J b , n a = n b , U a b/J = 12, V/J = 0.5, 
13 — L/J, and linear system L=12 (circles), 16 (squares), 
20 (up triangles), 24 (down triangles), 30 (diamonds). Both 
Sb (tv, 7r) and are non zero for a finite range of densities 

rib < 0.5, signifying the existence of a PSS phase. Errorbars 
are within symbol size when not visible. 

have found that, in order to observe PSF, one needs 
interactions U a b/J > 12. The relatively small effective 
hopping J c ff oc -jj— , or equivalently large effective 
mass, associated with PSF, and the presence of long 
range interaction renders the observation of a PSS a 
fairly challenging task. On one hand the superfluid 
stiffness is indeed suppressed by the enhancement of the 
effecticve mass of the composite object, on the other 
hand V/J s 1 gives rise to metastability [H)]. Since 
we have found that the CB solid remains robust against 
temperature and doping also close to the tip of the 
lobe, the optimal situation for the observation of PSS 
is indeed in this region of the phase diagram, where 
metastability is much less severe. 

In Fig. |4] we plot Sb{it,iT) (right axis) and n PSF (left 
axis) as a funciton of density with parameters J a = J b , 
n a = rif,, Uab/J = 12, V/J — 0.5, linear system sizes 
L=12 (circles), 16 (squares), 20 (up triangles), 24 (down 
triangles), 30 (diamonds), and inverse temperature 
j3 = L/J. There exist a finite range of densities for 
which both Sb{ir, n) and < (W a + Wb) 2 > are finite. 
(We have checked that S a (-K,n) and < (W a - W b ) 2 > 
are simultanuesly finite and zero respectively.) While for 
5b(7r, 7r) there are not finite size effects within errorbars 
and provided L > 12, for n p F finite size effects arc 
present. As L increases so does n PSF , signifying that 
in the thermodynamic limit n PSF is indeed finite. Pro- 
nounced size effects, which can be seen even at densities 
rif, ~ 0.43, where the CB solid is disappearing, are most 



likely due to low Kosterlitz-Thouless (KT) temperatures 
for the PSF-normal transition. This is further jutified 
by finite temperature results. We have indeed found 
that at n a = rib ~ 0.45, Tkt 0.05J (not shown here). 

The large effective mass of the composite object 
enhance metastability. As a result it is computationally 
challenging to establish how far the PSS phase extends. 
We have found the presence of many mctastable phases 
already at V/J ~ 1 and rib < 0.5, not very far from 
the tip. These states are characterized by the presence 
of various grain boundaries and/or extended defects. 
These observations are also relevant in the context of 
PSS with polar molecules in bilayered systems with no 
inter-layer tunneling [2i[. In [24| the authors study soft 
core bosons with in-plane nearest neighbor repulsion 
and inter-layer nearest neighbor attraction. By using 
a combination of mean field Gutzwillcr approach and 
exact diagonalization to study an effective hamiltonian 
of pairs, they predict the existence of a PSS phase. Polar 
molecules on bilayers with zero inter-layer tunneling 
can be interpreted as a two-component mixture with 
same intra-species long range interactions while the 
inter-species one, also long range, can be independently 
tuned (by changing the distance between layers). The 
model studied here differs from the one describing polar 
molecules in bilayers in that both, atom-atom and 
atom-molecule interactions are short range. We expect 
the metastability to become even more severe when 
all interactions are long range, as it happens in the 
case of bilayers, making the existence of a PSS phase 
questionable and surely extremely challenging to observe 
both in simulations and experimentally. Preliminary 
simulations validate this expectation. 

In conclusion, we have studied a mixture of hard core 
bosonic polar molecules and one atomic non reactive 
bosonic species, confined on a two-dimensional square 
lattice. We have shown that the presence of atoms 
renders the CB molecular solid phase more robust 
both at zero and finite temperature, i.e. CB forms at 
lower strengths of dipolar interaction and survives at 
higher temperatures. We have estimated that this effect 
could be sizeable in current experimental setups. As a 
byproduct, atoms also order in a solid phase with same 
melting temperatures as for molecules. We have also 
established that for large enough interaction between 
atoms and molecules a PSS can be stabilized. We 
expect the observation of the latter to be experimentally 
challenging due to the presence of metastability and the 
low temperatures at which the superfluid component 
disappears. 
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